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Abstract 

• i— I ' 

Supercontinuum generation by femtosecond filaments in air is investigated for different laser wavelengths ranging from ultraviolet 
Q to infrared. Particular attention is paid on the role of third-harmonic generation and temporal steepening effects, which enlarge the 
^ blue part of the spectrum. A unidirectional pulse propagation model and nonlinear evolution equations are numerically integrated 

and their results are compared. Apart from the choice of the central wavelength, we emphasize the importance of the saturation 

• ^ intensity reached by self-guided pulses, together with their temporal duration and propagation length as key players acting on 
£>^both supercontinuum generation of the pump wave and emergence of the third harmonics. Maximal broadening is observed for 

^ (-h large wavelengths and long filamentation ranges. 

Key words: Supercontinuum Generation, Femtosecond Filaments, Nonlinear Schrodinger Equation, Self-phase Modulation 
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1. Introduction 
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Third harmonic (TH) generation and supercontinuum 
C^) [SC) emission are two phenomena which have attracted 
broad interest in the past years [1,2,3,4,5]. An evident rea- 
O son is their direct application in atmospheric remote sens- 
" 55 m S measurements based on LIDAR (Light Detection And 
>~»Ranging) femtosecond laser setups [6]. In this context, spec- 
tral broadening originates from complex mechanisms that 
^^drive the long-range propagation of ultrashort pulses, when 
^ they form narrow filaments in optically-transparent media. 
The physics of isolated femtosecond filaments in air is 
nowadays rather well understood (see, e.g., [7] and refer- 
ee ences therein). It involves the competition between Kerr 
self-focusing and plasma defocusing, triggered whenever 
the input pulse power exceeds the critical power for self- 
focusing P cr ~ AQ/(27mo«2)- Here, Ao is the central laser 
wavelength, hq = 1 and ri2 are the linear and nonlinear re- 
fraction indices in air, respectively. For high enough powers, 
multiple filaments nucleated after an early stage of modula- 
tional instability have also been widely investigated [7,8,9]. 
They produce spectral patterns mostly analogous to those 
generated by a single filament, as filamentary cells emerge 
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in phase from the background field and possess the same 
phase link [10]. By comparing Terawatt (TW) multifila- 
mented beams with Gigawatt (GW) single filaments in air, 
this property was again verified in the UV-visible region 
(230-500 nm), where femtosecond self- focusing pulses cen- 
tered at 800 nm generically produce a tremendous plateau 
of wavelengths [11,12,13]. 

This latter phenomenon has recently become a subject of 
inspiration for several researchers. Two scenarios have been 
proposed for justifying the build-up of new wavelengths in 
the UV-visible range. On the one hand, temporal steep- 
ening phenomena undergone by the pump were shown to 
deeply modify the filament spectrum [14]. Full chromatic 
dispersion included in the optical field wave number k(uj) 
affects both the diffraction operator and the nonlinearities. 
This induces shock-like dynamics at the back edge of the 
pulse through space-time focusing and self-steepening ef- 
fects, which strongly "blueshift" the spectra. On the other 
hand, spectral broadening becomes enhanced by harmonic 
generation. The coupling of TH with an infrared (IR) pump 
produces a "two-colored" filament from pump intensities 
above 10 TW/cm 2 [15,16,17]. The amount of pump energy 
transferred into TH radiation depends on the linear wave 
vector mismatch parameter Afc = [3fc(cj) — fc(3w)] _1 fix- 
ing the coherence length L c = 7r/|Afc|. The smaller the co- 
herence length, the weaker TH fields. Along meter-range 
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distances, the TH component can stabilize the pump wave 
with about 0.5% conversion efficiency [13]. Experimental 
and numerical data reported ring structures embarking 
most of the TH energy and having a half-divergence an- 
gle of about 0.5 mrad [18]. This process contributes to cre- 
ate a continuous spectral band of UV-visible wavelengths 
[11,13,19]. 

Resembling spectral dynamics have also been reported 
from 1-mJ infrared pulses propagating in argon at atmo- 
spheric pressure, after subsequent compression by chirped 
mirrors [20]. Simulations of these experiments [21], discard- 
ing TH emission, revealed that temporal gradients inherent 
to the steepening operators are sufficient to amplify UV 
shifts and cover the TH bandwidth down to 250 and 210 nm 
for initial pulse durations of 10 and 6 fs, respectively. Very 
recently, numerical simulations [22] refound this tendency 
for atmospheric propagation, i.e., TH generation, while it 
affects the pump dynamics to some extent over long ranges, 
does not change significantly SC spectra, whose variations 
are mostly induced by the fundamental field in air. 

Despite these last results, we are still missing a detailed 
understanding of the key parameters which are supposed 
to drive SC generation. A first important parameter is, of 
course, the laser wavelength itself: How does the supercon- 
tinuum evolve when Ao is varied? This question was ad- 
dressed in Ref. [11] for various laser wavelengths, at which 
some spectral components were seen to merge. However, 
the model used a two-envelope approximation (for the 
pump and TH fields, separately). As emphasized in [22], 
splitting into TH and SC pump within envelopes becomes 
problematic when their respective spectra overlap inside 
a wide frequency interval where the basic validity condi- 
tion Aujj/ujj <C 1 (j = u),3w) may no longer be fulfilled. 
Actually, TH radiation produced through the nonlinear 
polarization needs to be described self-consistently from a 
single equation governing the total real optical field. This 
model was missing in Refs. [13,19], which made the role of 
TH overestimated in the white light emission. Another im- 
portant parameter is the length of the self-guiding range: 
Successive cycles of focusing and defocusing events pro- 
mote the creation of shorter peaks in the pulse temporal 
profile and lead to a maximal extension of the spectrum. 
A third potential player is the input pulse duration. In 
[21], this was shown to affect the spectra in noble gases 
for pulses containing a few optical cycles mainly. Clearing 
this aspect requires several simulations using distinct pulse 
durations and exploiting different propagation ranges. In 
connection, we demonstrate that spectral enlargements are 
directly linked to the level of maximum intensity: Steepen- 
ing operators as well as TH radiation broaden all the more 
the spectra as the intensity in the filament is high. 

The paper is organized as follows: Sec. 2 presents the 
model equations, namely, a unidirectional propagation 
equation for the total electric field that generates higher- 
order harmonics (mostly TH) through Kerr nonlincaritics. 
Results from this equation will be compared with those 
inferred from the "standard" nonlinear evolution equation 



(NEE) for the pump wave. The major difference between 
these two models lies in the production of the TH field and 
its coupling with the pump wave. Sec. 3 is devoted to the 
long-range propagation of 127- fs pulses in air described by 
the previous models. Emphasis is put on the influence of 
the central wavelength A (248, 800, 1550 nm). We discuss 
spectral modifications versus the height of / max , the input 
duration, together with the temporal steepening dynamics 
and merging between TH and pump spectral bands. Sec. 4 
revisits SC for short-range (focused) propagations. It is 
shown that / max becomes closer to analytical evaluations 
when the beam develops few focusing/defocusing cycles. 
In this configuration, a lesser broadening may be achieved. 
Sec. 5 finally summarizes the generic features resulting 
from our analysis. 

2. Models for pulse propagation and underlying 
physics 

Our unidirectional pulse propagation equation (UPPE) 
assumes scalar and radially-symmetric approximations. It 
also supposes negligible backscattering. These hypotheses 
hold as long as the beam keeps transverse extensions larger 
than the central laser wavelength and as the nonlinear re- 
sponses (together with their longitudinal variations) are 
small compared with the linear refraction index. Straight- 
forward manipulations of Maxwell equations allow us to 
establish the equation for the spectral amplitude of the op- 
tical electric field in the forward direction as [7] 

^=^VlE + ikH E + ^ NL , (1) 

where E(r,z,u) = (2n)^ 1 J E(r, z,t)e lUJt dt is the Fourier 
transform of the forward electric field component, z 
is the propagation variable, = r~ 1 d r rd r (r = 

x 2 + y 2 ) is the diffraction operator, ^ e = 1/c 2 , k(oj) = 
yl + {of)ijj/c is the wavenumber of the optical field 
depending on the linear susceptibility tensor x^ 1 K a de- 
fined at frequency u). In Eq. (1), .Fnl = Pnl + iJ/u is the 
Fourier transform of the nonlincaritics that include the 
nonlinear optical polarization Pnl and the current density 
J created by charged particles. Eq. (1) restores the earlier 
UPPE formulation proposed by Kolesik et al. [23] in the 
limit k\/k 2 {uj) < 1 {k\ = k 2 x + k$). 

For practical use, it is convenient to introduce the com- 
plex version of the electric field 

E = y/c[(S +£*), £ = -1= I ®{u)Ee- iojt duj, (2) 
v c i J 

where c\ = woMo/2fco employs the central wavenumber and 
frequency of the pump wave (ho = noWo/c) and 0(x) de- 
notes the Heaviside function. Because £ satisfies £*(ui) = 
£(—ll>)* (* means complex conjugate), it is then sufficient 
to treat the UPPE model (1) in the frequency domain uj > 
only. The field intensity can be defined by E 2 averaged 
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over an optical period at least. Expressed in W/cm 2 , it is 
simply given by the classical relation I = \£\ 2 . 

Concerning the nonlinearities, we assume a linearly 
polarized field. We consider a cubic susceptibility tensor 
keeping a constant value around ujq, so that Pnl 
contains the instantaneous cubic polarization expressed 
as P^(r,t) = eox^oE 3 . In addition, the phenomenon 
of Raman scattering comes into play when the laser field 
interacts with anisotropic molecules, in which vibrational 
and rotational states are excited. Depending on the tran- 
sition frequency f2i3 in three- level molecular systems and 
related dipole matrix element \i [24], the Raman response 
takes the form 

o (1) 2 rt , 
PRaman = ^r4r I ^~ SmC—^)E 2 (t>)dt> X E, 
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(3) 

where t\ — 1/u>r is the inverse of the fundamental rota- 
tional frequency and r 2 is the dipole dephasing time. Ex- 
pressed in terms of the rescaled complex field £ [Eq. (2)] 
and with appropriate normalizations [25] , Eq. (3) completes 
the cubic polarization as 
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(4b) 
(4c) 



where n 2 = 3x£ 3 V(4"o ce o) is the Kerr nonlinear index. 
Expression (4a) possesses both retarded and instantaneous 
components in the ratio xk- The instantaneous part ~ 5{t) 
of Eq. (4b) describes the response from the bound elec- 
trons. The retarded part <~ h(t) accounts for the Raman 
contribution, in which fast oscillations in E 2 give negligi- 
ble contributions, as n <~ t 2 ~ 70 fs are assumed to exceed 
the optical period <~ luq 1 . 

When free electrons are created, they induce a current 
density J = q e pv e , which depends on the electron charge 
q e , the electron density p and the electron velocity v e . J is 
computed from fluid equations involving external plasma 
sources and the electron collision frequency v e . At moderate 
intensities (< 10 15 W/cm 2 ), the current density obeys 
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d t J +VeJ =—E. 



(5) 



Assuming electrons born at rest, the growth of the electron 
density is only governed by external source terms, i.e., 



d t p = w(i)( Pnt 



(6) 



that include photo- ionization processes with rate W(I) and 
collisional ionization with cross-section 

*M= ,„ jl , ,. 2/ „ 2 v (7) 



TO e e n ci/ e (l + lu 2 /v 2 ) 

Here, p nt and C/j are the density of neutral species and the 
ionization potential, respectively. Electron recombination 



in gases is efficient over long (ns) time scales, and therefore 
we omit it. In Eq. (6), the rate for photo-ionization W(I) 
follows from the Perelomov, Popov and Terent'ev (PPT)'s 
theory [26] incorporating Ammosov, Dclone and Krainov 
(ADK) coefficients [27] (see also Ref. [28]). Optical field 
ionization theories stress two major limits bounded by the 
Keldysh parameter, 
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namely, the limit for Multi-Photon Ionization (MPI, 7 ^> 1) 
concerned with rather low intensities and the tunnel limit 
(7 < 1) concerned with high intensities, from which the 
Coulomb barrier becomes low enough to let the electron 
tunnel out. Here, E p denotes the peak optical amplitude. 
For laser intensities I = \£\ 2 < 10 13 - 10 14 W/cm 2 , MPI 
characterized by the limit 



W(I) -» Wmpi = o K I 



K 



(9) 



dominates, where K = mod([/i/fia>o) + 1 is the number 
of photons necessary to liberate one electron. The level of 
clamped intensity, I max , depends on the selected ionization 
rate. 

Energy lost by the pulse through single ionization pro- 
cesses is determined by a local version of the Poynting the- 
orem, yielding the loss current Ji OSSl such that Ji oss ■ E = 
UiW(I)(p nt — p). As a result, our UPPE model reads in 
Fourier space as 
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where 
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e lUJt dt. (11) 



In Eq. (10), Pnl and the expression containing the electron 
density p(r,t) [Eq. (6)] must be transformed to Fourier 
space, from which we retain only positive frequencies for 
the symmetry reasons given above. 

Alternatively, when a central frequency u>q is imposed, 
Eq. (1) restitutes the Nonlinear Envelope Equation (NEE), 
earlier derived by Brabec and Krausz [29]. We can make 
use of the Taylor expansion 
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k(uj) = kQ + k'uj + V, D = ^ 



k (n) 



n>2 



(12) 



wherew = u—uq, k' = dk/dLj\ u=aJo , = d n k/duj n \ UJ=UJo , 
and take the inverse Fourier transform of Eq. (1) in which 
terms with k(w) in their denominator are expanded up to 
first order in u) only. After introducing the complex-field 
representation £ = JJc lk ° z ^ w ° t , the new time variable 
t — > t — z jv g can be utilized to replace the pulse into the 
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frame moving with the group velocity v g = k' 1 . Further- 
more assuming v^jj 1 <C 1, y/e(uja) / e(u)) w 1 and ignoring 
the TH component, the nonlinear envelope equation for 
the forward pump envelope U expands as follows: 

d 



^-U = ^-T-^lU + WU + i—n 2 Tx 
oz 2k c 

(l-x K )\U\ 2 +x K [ h(t-t')\U(t') 

J —oo 



dt' 



u 



2\U\< 



(13) 



where a = a(u ), V = ES( fc( "V« ! )(^t)" an d T = 
(1 + ^d t ). The first term of the operator T> corre- 
sponds to group-velocity dispersion with coefficient k" = 
d 2 k/duLi 2 \ lx j =U j Q . Equation (13) describes wave diffraction, 
Kerr focusing response, plasma generation, chromatic dis- 
persion with self-consistent deviations from the classical 
slowly- varying envelope approximation through the space- 
time focusing and self-steepening operators [(T _1 Vj_£) 
and (T\£\ 2 £ ), respectively]. This model will be integrated 
numerically by using initially Gaussian pulses, 



U(x,y,z = 0,t) 
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which involves the input power P m , the beam waist wq 
and 1/e 2 pulse half- width t p . Input pulses can be focused 
through a lens of focal length / and they linearly diffract 
over the distance 

z f = {f/z )/{l + f/zl), (15) 

where zq = Tmow 2 /\o is the Rayleigh range of the colli- 
mated beam (/ = +oo). 

In the coming analysis, wc shall employ the nonlin- 
ear refractive indices n 2 = 8 x 10~ 19 , 4 x 10~ 19 , and 
1 x 1CT 19 cm 2 /W for the wavelengths A = 248, 800 and 
1550 nm, respectively. At 800 nm, we consider a fitted 
MPI formulation for the ionization rate, W(I) — > <j/k)I K j 



where K = 8 and cr (8) = 2.88 x 10"" s^cm^/W 8 . This 
approximation is known to reproduce experimental data 
at A = 800 nm rather faithfully [28,30]. For the two other 
wavelengths, we lack well established formulations, so we 
employ PPT ionization rates. We consider 2 molecules 
having the lowest gap potential (C/j = 12.1 eV) as the 
main specy undergoing ionization with an effective resid- 
ual charge Z e g = 0.53 [31]. All ionization rates used in 
the present paper are illustrated in Fig. 1. They yield 
saturation intensity below the threshold of 100 TW/cm 2 
currently claimed in the literature [7,11]. The dispersion 
relation for air has been parametrized as in Ref. [32] . 

Since our outlook is to understand spectral variations 
versus the propagation dynamics at different wavelengths, 
we find it instructive to fix the same ratio of input power 
over critical, e. g., P; n = 4x P cr , at all wavelengths. To lo- 
cate the Kerr-driven filamentation onset upon comparable 
z scales, we also adjust the ratio of the nonlinear focus z c 
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Fig. 1. Ionization rates for the different wavelengths used throughout 
this paper: Ao = 800 nm (solid line), Ao = 248 nm (dotted line), 
and Ao = 1550 nm (dashed line). The dashed-dotted line shows the 
overestimated ionization rate chosen in Sec. 3.2 for Ao = 1550 nm. 

and Rayleigh range z between 2 and 4, by adapting suit- 
ably the beam waist wo between 1 and 4 mm. 

2.1. Self-phase modulation and SC generation 

In air, dispersion is weak with k" < 1 fs 2 /cm in the 
UV as well as in the mid-IR domains. Self-channeling then 
mainly relies on the dynamical balance between Kerr self- 
focusing and plasma defocusing, so that estimates for peak 
intensities (I max ), electron densities (/0 max ) and filament 
radius (L m ; n ) can be deduced from equating diffraction, 
Kerr and ionization responses in Eq. (13). This yields the 
simple relations 
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2p c n n 2 

L m in ~ 7r(2fc 2 n 2 /ma X /«o) _1/2 , 



(16a) 
(16b) 



where 



"2 = "2(1 — xk) + n 2 a;xmaxi / h(t — t')e '? dt' ', (17) 

J — 00 

represents the maximal effective Kerr index over the initial 
pulse profile. For practical use, W(I max ) can be simplified 
to aKl^ ax in MPI-like formulation. 

The magnitude of J max directly impacts spectral broad- 
ening, which is basically driven by self-phase modulation 
(SPM). Because the frequency spectrum is expanded by 
the nonlinearity, SPM leads to SC, as the wave intensity 
strongly increases through the self-focusing process. Noting 
by <p(r, t) the phase of the field envelope, frequency varia- 
tions are dictated in the limit T, T^ 1 — > 1 by 

Auj = -d t <p ~ -k Azdt(n 2 I - p/2n p c ), (18) 

which varies with the superimposed actions of the Kerr and 
plasma responses. Near the focus point z c , only the front 
edge of the pulse survives from this interplay and a red- 
shift is enhanced by plasma generation. At later distances, 
second focusing/defocusing sequences attenuate this first 
tendency. In contrast, when accounting for temporal steep- 
ening (T, T _1 7^ 1), shock edges in the back of the pulse are 
created and a "blue shoulder" appears in the spectrum, to 
the detriment of the early redshift [7,14,33,34]. 
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In addition, the cubic polarization generates third-order 
harmonics, modeled by the last term of Eq. (4a). In self- 
focusing regimes, the third-harmonic intensity usually con- 
tributes by a little percentage to the overall beam fluence 
[16]. Despite its smallness, this component may act as a sat- 
urable nonlinearity for the carrier wave. It lowers the peak 
intensity of the pump and contributes to enhance the blue 
side of the spectrum after the TH and pump bandwidths 
increase and overlap [11,13]. 



3. Long-range Propagation 

We numerically analyze supercontinuum generation for 
the three laser wavelengths of 248 nm, 800 nm and 1550 
nm. Results of the nonlinear Schrodinger-like equation (13) 
for the pump envelope involving or not space-time focusing 
and self-steepening operators are compared with those of 
the unidirectional propagation equation (10) avoiding any 
Taylor expansion in the dispersion relation. 

Special attention is here given to the long-range prop- 
agation, which rather favors several cycles of focusing- 
defocusing events. Before proceeding with the above pa- 
rameters specifically, we perform three different series of 
simulations showing SC at 800 nm, whose results are sum- 
marized in Fig. 2. The first one concerns direct integrations 
of Eq. (10); the second refers to the same pulse described 
by Eq. (13), which eludes TH production; the third ap- 
proach relies on Eq. (13), in which temporal steepening is 
omitted, i.e., T = T~ 1 = 1. 

The insets in Fig. 2 show the on-axis spectra when SC 
is maximal. Following the UPPE description, TH, which 
emerges from z c ~ 6 m, develops a limited redshift, whereas 
SC of the fundamental is widely extending towards the 
blue/UV wavelengths [Fig. 2(a)]. Note that, although a 
broad plateau occurs in this domain, the TH bandwidth 
still appears separated from the pump spectrum. Following 
the NEE description, there is no TH generation. However, 
SC is so amplified in the blue region by temporal steepen- 
ing effects, that it overlaps the TH zone and simply hides 
it [Fig. 2(b)]. Finally, when neglecting temporal steepen- 
ing, the pump instead develops a wide redshift (overesti- 
mated by plasma coupling) and a much narrower blucshift 
[Fig. 2(c)]. 

A first observation can be drawn from Fig. 2: Since TH 
is responsible for lowering the saturation intensity of the 
pump [13], / max reached in the UPPE model is lower and 
the frequency variations Alo <~ J max Az/Af are diminished 
compared with NEE spectra for the pump wave alone. 
Apart from this difference, no significant other change was 
detected between both these models, so that NEE seems 
to be nothing else but the UPPE description subtracted 
by the self-generated harmonics [7] . Importantly, omitting 
temporal derivatives of the operators T, T _1 imply more 
serious discrepancies, as can be seen from Fig. 2. 
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Fig. 2. Peak intensities (solid curves, left-hand side scales) and peak 
electron densities (dashed curves, right-hand side scales) of 127-fs 
pulses with ratio of input power over critical equal to 4 and waist 
wo = 2 mm at Aq = 800 nm. The insets show on-axis spectra: dotted 
curves represent the initial spectra; the solid curves the spectra at 
the propagation distance z m ax = 8 m where maximal broadening is 
observed: (a) UPPE model Eq. (10); (b) NEE Eq. (13) applied to the 
pump wave; (c) NEE Eq. (13) modified with setting T = T" 1 = 1. 

3.1. Influence of Ao 

We now examine Gaussian pulses at different wave- 
lengths (248, 800 and 1550 nm) with t p = 127 fs and 
Pin/ Per = 4 as initial conditions for the UPPE model in 
parallel geometry (/ = +oo). Figures 2(a) and 3 show 
snapshots of spectra at maximal extent, together with as- 
sociated peak intensities and electron densities. We here 
specify that no TH generation was included for A = 248 
nm, because no reliable data of the dispersion relation 
was available for this wavelength. We believe, instead, 
that spectral components below 90 nm should be rapidly 
absorbed by the medium. 

By comparing Figs. 2(a) and 3, it is seen right away that 
supercontinuum generation increases with the wavelength. 
To quantify this observation we introduce AAsc as the to- 
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Fig. 3. Peak intensities (solid curves, left-hand side scales) and peak 
electron densities (dashed curves, right-hand side scales) of 127-fs 
pulses with ratio of input power over critical equal to 4 at differ- 
ent wavelengths Ao- The insets show on-axis spectra: dotted curves 
represent the initial spectra; the solid curves the spectra at the 
propagation distance z ma x where maximal broadening is observed: 
(a) Ao = 248 nm, wo = 1 mm, z ma x = 6 m; (b) Ao = 1550 nm, 
wo = 4 mm, z max = 35 m. 
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tal extension of the on-axis spectra over wavelengths at 
10 -5 times the maximal spectral intensity. Then a mea- 
surement for the effective broadening is the ratio AAsc / -V), 
which we find close to ~ 0.5 at 248 nm, ~ 1 at 800 nm 
and ~ 1.5 at 1550 nm. A look at the propagation dynam- 
ics reveals that the self-guiding range is noticeably aug- 
mented at longer wavelengths. This can be explained by 
the transverse size of the filament. Equation (16b) gives 
an estimate for the beam waist in filamentation regime. 
If we assume comparable I max for all wavelengths, we de- 
duce that the filament diameter at 1550 nm is about one 
order of magnitude larger than that at 248 nm, which is 
compatible with our numerical data. Hence, the larger the 
wavelength, the slower the filament is expected to diffract. 
Moreover, by virtue of the formula for the critical power 
P CT oc Ag/ri2 and since a filament conveys a few P cr [7], it is 
obvious that IR filaments contain much more energy than 
their UV counterparts. Thus, nonlinear losses along the fil- 
amentation range are less dramatic in the IR domain. Due 
to the longer propagation range, more focusing/defocusing 
events participate in enlarging the spectra at longer wave- 
lengths. Figure 4 details the evolution of the filaments in 
the plane (t, z). It is seen that the time window in which the 
pulses disperse occupies the length of the input pulse dura- 
tion. Although shorter temporal peaks arise through self- 
focusing/defocusing events, multi-peaked profiles mostly 
develop patterns having a whole extent close to t p . 
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Fig. 4. Temporal evolutions of the pulses shown in Figs. 2(a) and 3 in 
the (t, z) plane: (a) A = 248 nm; (b) A = 800 nm; (c) A = 1550 nm. 

Before going on, we find it worth investigating SC at 1550 
nm more thoroughly, in relationship with third-harmonic 
generation. Figure 5 plots three of the SC development 
stages, first when TH and pump components are clearly 
separated (z = 10 m), second when they start to merge 
(z = 15 m). At the last propagation distance [z = 25 m), 
we can observe that, unlike Fig. 2, TH and pump spectra 
overlap and make the TH bandwidth not distinguishable. 
For comparison, results from the NEE model for the pump 
wave alone have also been plotted. 

3.2. Influence o// max 

The impact of the saturation intensity onto SC is in- 
vestigated by simply changing the ionization model: By 
decreasing the photo-ionization rate artificially it is pos- 
sible to increase / max and the maximal plasma level p m ax 
accordingly. Reversely, increasing W(I) reduces these two 
quantities, which can have a direct influence on the spec- 
tral broadening, as inferred from Eq. (18). To study this 
point, we concentrate on the wavelength of 1550 nm only, 



G 




10" 



_ 10"' 

<< 

o" 

n 

- 10" 4 



10" 



0.5 


1 1.5 2 2.5 

X (itm) 


(b) 






















IPwl ' 
















1 i 






1 i 







0.5 1 



1.5 2 

X (itm) 



2.5 



Fig. 5. On-axis spectra for the pulse used in Fig. 3(b) (Ao = 1550 nm) 
for the propagation distances (a) z = 10 m (solid line) and z = 15 m 
(dash-dotted line) and (b) z = 25 m (solid line). The dashed line in 
(b) refers to a spectrum at z = 25 m computed from the NEE model. 

because it yields the broadest spectra explored so far. Fig- 
ure 6 shows the maximal intensity and peak electron den- 
sity at this wavelength, when the ionization rate is artifi- 
cially increased (see Fig. 1). All other parameters are un- 
changed, compared with the simulation shown in Fig. 3(b). 
With the original ionization rate, 7 max reaches the value 
of 80 TW/cm 2 ; with the artificial one, 7 max stays below 
13 TW/cm 2 . The inset in Fig. 6 details the corresponding 
spectrum at z max = 40 m. With low I max , the TH compo- 
nent is reduced to some extent as the pump intensity barely 
exceeds the TH conversion threshold [15]. Meanwhile, SC 
of the pump driven by the T, T _1 operators in the blue side 
decreases, i.e., a lower / max for analogous pulse compres- 
sion implies less sharp temporal gradients and smoother 
optical shocks, which weakens blucshifted frequency vari- 
ations. These features are visible in Fig. 6, where TH and 
pump broadbands remain separated. In Figs. 3(b) and 5(b), 
in contrast, SC extends beyond the TH wavelength and in- 
creases more the pulse spectrum. Similar features were ob- 
served at the two other wavelengths, when the ionization 
rate was changed. 



3.3. Influence oft p 

We investigate the influence of the initial pulse duration 
on the propagation dynamics and SC generation. Since we 
consider transform- limited pulses, the value of t p is directly 
linked to the initial spectral width. Moreover, J max comes 
into play in SC generation and is expected to scale as ~ 
(l/fp) 1 /^" 1 ) [sec Eq. (16a)]. Thus, the initial pulse dura- 




40 60 
z(m) 



Fig. 6. Peak intensity (solid curve, left-hand side scale) and peak 
electron density (dashed curve, right-hand side scale) for the same 
pulse as in Fig. 3(b) (Ao = 1550 nm) computed from the UPPE 
model with an overestimated ionization rate (see Fig. 1). The inset 
shows maximal spectral broadening attained at z max = 40 m. 

tion should play a significant role in spectral broadening. 
To check this assessment, we performed several simulations 
using the UPPE model, by varying t p from 20 fs up to 
500 fs. Because group- velocity dispersion becomes very ef- 
ficient at short pulse durations and may even stop the Kerr 
self-focusing at powers too close to critical [35,36], we in- 
creased the input power up to 20 P CI for t p = 20 fs. With 
this, we ensure to trigger a filamentation regime even for 
this short input duration. 

To illustrate the dependency of I max upon t p , short wave- 
lengths are preferable because the number of photons for 
ionization is small. Figure 7 shows maximum intensity, peak 
electron density and maximal spectral extent of a 20-fs 
pulse at 248 nm. At this wavelength, K = 3 and, follow- 
ing Eq. (16a), 7 max and p max should increase by a factor 
<~ 2.5 compared to the 127-pulse shown in Fig. 3(a). In- 
deed, both quantities are increased by a factor of ~ 2 in 
the simulation. Spectral broadening is augmented from 0.5 
to 0.8 in terms of AAsc/Ao, especially to the blue side. 
As explained in Sec. 3.2, this results from the action of 
the steepening operators. The overall propagation dynam- 
ics, characterized by the filamentation length and number 
of focusing/defocusing cycles are, however, comparable for 
both the 20-fs and the 127-fs pulses. 

On the other hand, if we increase the initial duration 
t p towards the ps time scale, the propagation dynamics 
changes drastically. As an example, Figure 8 shows the 
temporal evolution of a 500-fs pulse at 800 nm. Compared 
with Fig. 4(b) employing t p = 127 fs, the obvious difference 
is the huge number of focusing/defocusing cycles. The ac- 
tion of the generated plasma breaks the pulse profile into 
a larger number of shorter peaks. With a longer pulse du- 
ration, more "time slices" are available for feeding succes- 
sive focusing events. The filamentation range is increased 
and 7 max is maintained over several meters. Inspection of 
the simulations, however, reveals maximal spectral extent 
comparable with that dispayed in Fig. 2(a). 

At 1550 nm the influence of t p on the maximal inten- 
sity is much less pronounced, since we have K = 15 [see 
Figs. 3(b) and 9]. For all pulse durations, we indeed ob- 
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Fig. 7. Peak intensity (solid curve, left-hand side scale) and peak 
electron density (dashed curve, right-hand side scale) of a 20-fs pulse 
with ratio of input power over critical equal to 20, Aq = 248 nm, 
wo = 1 mm. The inset shows on-axis spectra: the dotted curve 
represents the initial spectrum; the solid curve the spectrum at the 
distance z m ax = 3.5 m where maximal broadening is observed. 
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Fig. 8. Temporal evolution of a 500-fs pulse at 800 nm with 4 P c , 
and waist wo = 2 mm in the (t, z) plane. 



serve I max ~ 80 TW/cm 2 . This can explain why the max- 
imal spectral extension AAsc/^o is always found between 
1.5 and 2, regardless t p may be. The major difference lies 
in the filamentation range, which increases with the initial 
pulse duration. So, there is no significant change in SC gen- 
eration between short and long pulses over large enough 
propagation scales. 

This last observation invites us to look at the temporal 
profiles upon propagation. If the spectral extent is compa- 
rable, we might also find similar temporal patterns. Indeed, 
the on-axis temporal profiles shown in Fig. 10 all exhibit 
structures with duration of 10-15 fs. It seems that the initial 
pulse length t p just determines how many of these peaks 
appear, or, in other words, how many focusing/defocusing 
cycles the pulse is able to develop upon propagation. An- 
other indication for the change in the effective pulse dura- 
tion upon propagation is provided by the curve of the in- 
tensity maximum in Fig. 9(b): The first focusing cycle is 
halted at slightly lower intensities ~ 60 TW/cm 2 , because 
in this early stage the peak duration remains of the order 
of t p = 500 fs. At later stages, J max increases as the pulse 
undergoes temporal compression. 




Fig. 9. Peak intensities (solid curves, left-hand side scales) and 
peak electron densities (dashed curves, right-hand side scales) of 
4 mm waisted pulses with different durations and ratios Pin/Per at 
Ao = 1550 nm. The insets show on-axis spectra: dotted curves rep- 
resent the initial spectra; the solid curves the spectra at the prop- 
agation distance z m ax where maximal broadening is observed: (a) 
t p = 20 fs, P in = 15 X P CI ., z max = 20 m; (b) t p = 500 fs, P in = 4xP cr , 

Zmax = 40 m. 
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Fig. 10. On-axis temporal profiles of pulses with different duration at 



: 500 fs [dotted line, 



[dashed line, parameters used in Fig. 3(b)]; t p 
parameters used in Fig. 9(b)]. 

4. Short-range Propagation 



So far, we have analyzed free propagation dynamics 
where long filaments achieve temporal gradients and SC 
extents similar to those produced by initially much shorter 
pulses. Now, we force all pulses to cover the same short 
filamentation range through a focusing optics (/ = 2 m). 
Results are shown in Fig. 11 for t p — 20 and 500 fs at 
1550 nm. We can check that the maximum intensity J ma x 
follows the theoretical expectations (16) involving t p . The 
short filamentation range prevents the occurrence of several 
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Fig. 11. Peak intensities (solid curves, left-hand side scales) and 
peak electron densities (dashed curves, right-hand side scales) of 
4 mm waisted pulses with different duration and ratio of input power 
over critical equal to 4 at Ao = 1550 nm propagating in focused 
geometry (/ = 2 m). The insets show on-axis spectra: dotted curves 
represent the initial spectra; the solid curves the spectra at the 
distance z ma x = 2 m where maximal broadening is observed: (a) 
t p = 20 fs; (b) t p = 500 fs. (c) and (d) show the respective temporal 
evolution for the pulses of (a) and (b) in the (t, z) plane. 



focusing/dcfocusing cycles, especially for the 500-fs pulse. 
Hence, although the pulse self-focuses, its temporal extent 
remains comparable with t p , as we can see in Figs. 11(c) 
and (d). The 20-fs pulse shows significant spectral broad- 
ening with a visible, broadened TH peak. The 500-fs pulse 
is spectrally too narrow to generate supercontinuum over 
roughly 1-m filamentation range. Therefore, the funda- 
mental and harmonic peaks clearly stay separated (note 
the occurrence of the fifth harmonics at 310 nm, which is 
self-consistently described by the UPPE model). Thus, the 
initial pulse duration strongly influences spectral broaden- 
ing in configurations of short filamentation range mainly, 
which is consistent with the numerical results of Ref. [21]. 



5. Conclusion 

In summary, we have revisited recent works on SC gen- 
eration versus third harmonic emission, by showing from a 
complete UPPE model that spectral enlargements of fem- 
tosecond pulses in self-guiding regime are mostly driven by 
space-time focusing and self-steepening. TH generation, al- 
though changing the pump dynamics, affects the spectra to 
a limited extent only. This conclusion corrects some previ- 
ous statements [19] , based on a propagation model in which 
temporal steepening terms were analyzed separately from 
an envelope description for TH generation. Going one step 
beyond, we have demonstrated the important role of the 
saturation intensity in the frequency variations enlarging 
both TH and pump broadbands. The input pulse duration 
becomes a significant player in the spectral extents as long 
as pulses do not propagate too far, i.e., they do not let the 
temporal profiles of the pulse fluctuate so much that many 
peaks and sharp gradients can develop along the optical 
path. This property is mainly enlightened in short-range 
focused geometry and lost in long-range parallel geome- 
try. Finally, the role of the central wavelength is preem- 
inent: Our numerical simulations displayed evidence that 
SC clearly augments with the laser wavelength. 
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